Global strong solution to a nonlinear Dirac type 
equation in one dimension 

Yongqian Zhang 

School of Mathematical Sciences, Fudan University, Shanghai 200433 
Key Laboratory of Mathematics for Nonlinear Sciences 
Email address: yongqianz@fudan.edu.cn 

March 4, 2013 



Abstract 

This paper studies a class of nonlinear massless Dirac equations in 
one dimension, which include the equations for massless Thirring model 
and massless Gross-Neveu model. Under the assumptions that the ini- 
tial data has small charge and is bounded, the global existence of the 
strong solution is established. The decay of the local charge is also 
proved. 

1 Introduction 



We consider the nonlinear massless Dirac equations 

f i(u t + u x ) = Gi(u,v), 
1 i(v t - v x ) = G 2 {u,v), 



[1.1) 



with initial data 

(u,v)\ t =o = {u {x),v {x)). (1.2) 

where (x, t) G R 2 , (u, v) G C 2 , Gq = daW(u, v)+F 1 (u, v) and G 2 = ckW{u, v)+ 
F2(u,v) which satisfy the following: 

(Al) W(u,v) : C 2 i— > R is a polynomial in variables \u\ 2 and \v | 2 ; 

(A2) Fk (k = 1,2) are polynomials of cubic order in variables u,v and u,v, 
moreover there exists a positive constant c such that 



\uFi(u,v)\ + \vF2(u,v)\ < c\u\ 2 \v I 2 



for (u,v) G C 2 . 
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Many physical models verify (Al) and (A2), for examples, 

• Massless Thirrring model ([22]) 

W = a\u\ 2 \v\ 2 , a G R, F 1 = F 2 = 0; 

• Federbusch model 

W = a(\u\ 2 + |u| 2 )|u| 2 |u| 2 , a G R, F l = F 2 = 0; 

• Massless Gross- Neveu model ([15]) 

(Fi, F 2 ) = W { u^a(uv + uv) 2 , a G R, W = 0. 

The nonlinear Dirac equation has been studied by many authors for the mass- 
less or massive models, see for examples [1], [2], [5], [6], [9], [12], [10], [11], [13], 
[17], [18], [19], [20], [22] and the references therein. The global existence of so- 
lution for Cauchy problem has been established in [1], [9], [10], [11] in Sobolev 
space H s for suitable s. Recently some studies are devoted to the local low 
regularity solutions in a certain subspace of C([0, T]; H S (R)) for s > 0; see 
for instance [18], [19] and [20] and references therein. In [17], using the explicit 
representation of the smooth solution by [11] (see also [18]), Huh got the strong 
solution in L 2 (R) of the nonlinear equations for the massless Thirrring model 
and Federbusch model. The uniqueness in L°°([0,T] : L 2 (R) flL 4 (i?)) and the 
asymptotic behavior of the field were also established there. In [5], Candy used 
the strichatz estimates to prove global existence of L 2 -strong solution to the 
nonlinear Dirac equation for massive Thirring model. For the nonlinear term 
being a polynomial of quintic or more higher order, the small global solution 
has been constructed in [19] by the end point Strichartz estimates. 

Different from the above work, we consider in this paper the case that 
nonlinear terms include some cubic terms and may have the more complex 
structure. Such kind of equations include the Massless Gross-Neveu model. To 
our knowledge, the global existence of strong solution in L 2 for Gross-Neveu 
model is still open [19]. Different from the Thirring model, we have (2.1) for \u\ 2 
and \v\ 2 instead of the conservation law for the charge /^(M 2 + \v \ 2 )dx. The 
assumption (A2) implies that the equation (2.1) is similar to the approximate 
conservation law for the interaction between the elementary waves [14]; see 
also [4] [7] and [21]. Following Glimm's approach, we introduce a Glimm type 
functional to get the L 2 and L°° bounds of the local H 1 solution for extending 
the solution globally in time. Such functional was first given by Glimm [14] to 
establish the global existence of the BV solution for the system of conservation 
laws, and similar technique has been used to study the general conservation 
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laws, Boltzmann equations and wave maps in R 1+1 , see for instance [4], [7], 
[16], [14], [21] and [23] and the references therein. Our aim is to find the global 
solutions in the following sense. 

Definition 1.1 A function (u,v) is said to be a strong solution to (1.1-1.2) on 

the time interval [0,T] if there exists a sequence of classical solutions (u n ,v n ) 
to (1.1) on [0,T] x R such that (u n ,v n ) ->■ (u,v) strongly in C([0, T];L 2 (R)) 
and (u n ,v n )\ t=0 — > (u ,v ) strongly in L 2 (R). 

The main results is stated as follows. 

Theorem 1.1 Suppose that (uq,vq) G L 2 (i? x ) D L 00 ^ 1 ) and that the norm 
I |i*o I \l 2 (r) + H^olU 2 is small. Then (1.1-1.2) admits a unique global strong 
solution (u,v) G C([0, oo); L 2 (i? 1 )) f] L°°(R l x [0, oo)). Moreover, for any 
constants a <b, 



The remaining part of the paper is organized as follows. In section 2, based 
on the existence of local if 1 — solutions, we will introduce a Glimm functional 
for H l solutions and get the uniform L 2 and L°° bounds of the solution, which 
are independent of the time and yield the bound of H 1 norm for the solution. 
This enable us to extend the smooth solutions globally in time; see Proposition 
2.1. In section 3, we take the limit to get the global strong solution in L 2 and 
prove the decay of the local charge. 

2 Global classical solution 

We consider the case that (u ,v ) G in this section. (1.1) is a semilinear 
hyperbolic system. Local existence of the solutions in H 1 can be proved by the 
standard methods using the Duhamel's principle and the fixed point argument, 
see for instance [2], [9] and [21]. Now suppose that for (u , v ) £ C£°> the initial 
value problem (1.1-1.2) has a smooth solution (u,v) G C([0, T*); if 1 (i? 1 )) fl 
C 1 ([0, T*); L 2 (R)) for some T* > 0. Then the support of (u(t,x),v(t,x)) is 
compact for each t G [0, T*). To extend the solution across the time t = T* we 
only need to show that 



where the constants C and C" are independent of T*. To this end, we multiply 
the first equation of (1.1) by u and the second equation by v, then 




sup (Ht)|| ff i + \\v(t)\\ H i) < C'exp(C"n) 



0<t<T» 



(\u\ 2 ) t + (\u\ 2 ) x = 2^F lU ), 
(\v\ 2 ) t -(\v\ 2 ) x = 2^F 2 v). 



(2.1) 
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Here and in sequel tflZ denotes the real part of Z e C. 

As in [3] (see also [16] and [23]), we define the following functional for the 
solution for any t G [0,T*), 



lx<y 

and 



Q(t) =11 \u{x,t)\ 2 \v{y,t)\ 2 dxdy, 

J J x<y 



L{t) = / (\u(x,t)\ z + \v(x,t)\*)dx, 

J — oo 

/oo 
\u(x,t)\ 2 \v(x,t)\ 2 dx. 
-oo 

These functionals were given in [3] (see also [16]). Q(t) is called Bony 
functional [16] and is similar to the Glimm interaction potential. The lin- 
ear combination L(t) + KQ(t) for any constant K > is called the Glimm 
functional [14]. Direct computation shows that 



dQjt) 
dt 



< 



< 



-ff {\u{x,t)\ 2 ) x \v{y,t)\ 2 dxdy+ [ [ \u{x,t)\ 2 {\v{y,t)\ 2 ) y dxdy 

J J x<y J J x<y 

+2 / / \uF 1 (u,v)(x,t)\\v(y,t)\ 2 dxdy 
+2 f j \u(x,t)\ 2 \vF 2 (u,v)(y,t)\dxdy 

J J x<y 

/OO 
t)\ 2 \v(x,t)\ 2 dx 
-oo 

/oo roo poo poo 

\uFi(u,v)\dx / \v\ 2 dy + 2 / |w| 2 <ir / \vF 2 (u,v)\dy 
-oo J— oo J —oo J —oo 

/oo 
\u(x,t)\ 2 \v(x,t)\ 2 dx 
-oo 



Therefore for any constant X > 0, 

+ ^^(0 ^ {( _ 2 + 2L(t)c) ^ + 2c]D{t y (2 2) 
at at 

Lemma 2.1 There exzsts constants 5 > and > s«c/i t/iat /or t/ie initial 
data satisfying L(0) + fT<5(0) < 5 t/iere ZioZds t/ie following 

dL(t) dQ{t) , . 

-sr + * -Tb^m (2 ' 3) 



and 
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fort G [0, T*). Therefore, 

L{t) + #Q(i) + / T £>(r)dr < L(0) + iTQ(O) (2.4) 
Jo 

forte [0,T,). 

Proof. First choose constants 5 and X so that 

-2 + 25c < -1/2, -fT/2 + 2c<-2. 

If L(0) + XQ(0) < 5, then (-2 + 2L(0)c)ir + 2c < -2. Let 

T = sup{T; (-2 + 2L(t)c)K + 2c < —1, t G [0, T)}. 

We want to show that T = T*. Otherwise, suppose that T < T*, then with 
(2.2) we have (2.4) for t G [0,T ). Therefore, by the continuity of L(t) and 
Q(t) for < t < T*, we have the following 

L(t) + KQ(t) < 5, (-2 + 2L(t)c)K + 2c< -2, ie[0,T ], 

which leads to the contradiction to the definition of T by using again the 
continuity of L(t) for < t < T*. Thus the proof is complete. □ 

Next we shall use the above estimates to get the L°° bound of the solution 
for < t < T*. Denote 

n(x , t ) = {{x, t) : < t < t , x - (t - t) < x < x + (t - t)}, 

and 

T R (x , x Q ) = {(x, t) : < t < t , x = x + (t - t)}, 
r L (x , x ) = {(x, t) : < t < t , x = x - (t - t)}, 

see Fig. 1. 

Lemma 2.2 For any t G (0, T*), there hold 

r — r r i rxo+to—t 

/ \u\ 2 <V2c / \u\ 2 \v\ 2 dxdt + —= \ |m (^)| 2 ^ 

Jrri(xo,t ) J JQ(xo,to) \/2 Jxo—to+t 

and 

f f f | fXO+tQ— t 

/ \v\ 2 < V2c / / \u\ 2 \v\ 2 dxdt + —= \ \vo(x)\ 2 dx. 

JT L (x ,t ) J Jn(x ,t ) v2 Jx -t +t 
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0&o,*o) 




Figure 1: Domain Q(xo,to) 



Proof. Integrating (2.1) over Q(xo,to), then by Green formula we have 

rxo+to-t 



V2 \u\ 2 = / / 2$t{iF lU )dxdt + / \u (x)\ 2 dx 

r r rxo+tQ—t 

< 2c \u\ 2 \v\ 2 dxdt + \u (x)\ 2 dx 

J JQ(xci,ta) Jxa—tn+t 



and 

rxo+to-t 



V2 f \v\ 2 = f f 2^{iF 2 v)dxdt+ f X ° ° \v {x)\ 2 dx 

Jr L (xo,to) J JSl{xQ,to) Jxo-to+t 

r r rxo+to—t 

< 2c \u\ 2 \v\ 2 dxdt + / \v (x)\ 2 dx. 

J J(l(xa,ta) Jxn-tn+t 



The proof is complete. □ 



Lemma 2.3 Let 5 and K be the constants given in Lemma 2.1. Then for any 
constants a and b with a < b there hold the following 

sup (\u(x,t)\ 2 + \v(x,t)\ 2 ) <d sup (\u (x)\ 2 + \v (x)\ 2 ), 

0<t<T t ,a<x<b a-T*<x<b+T* 

and 

f (\u(x,t)\ 2 + \v(x,t)\ 2 )dx < Ci f (\u (x -t)\ 2 + \v (x + t)\ 2 )dx, 

J a J a 

where d = exp((2 v / 2c 2 + y/2c)S). 
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\u 



(x + (t-t ),t)\ 2 = 2%t(iF lU ) 



Proof. Using the characteristic method, by (2.1) and (A2) we have 
d 

< 2c\u(x + (t-t ),t)\ 2 \v(x + (t-to),t)\' 

and 

d_ 

dt' 

< 2c\u(xo-(t-to),t)\ 2 \v(x -(t-t ),t)f 
for any (x , t ) with a — (T* — t ) < x < b + (T* — t ). Then, 

rto 

\u(x ,t )\ 2 < \u(x - t ,0)| 2 exp (2c / \v(x + (r - t ), r)\ 2 dr) 



\v(x -(t-t ),t)\ 2 = 2M(iF 2 v) 



and 



= Kx -t ,0)| 2 exp(2c ( \v\ 2 ) 

Jr L (x ,t ) 

< \u (x - t )| 2 exp (h(x ,t )) 



rto 

\v{x ,t )\ 2 < \v(x + t ,0)\ 2 exp(2c \u(x - (r - t ),T)\ 2 dr) 

Jo 

= \v(x + t ,0)\ 2 exp(2c f \u\ 2 ) 
< |v(xo + to,0)| 2 exp(/ 2 (a;o,to)), 



where the estimates for Jr R ( Xo ,t ) \ u \ 2 an d fr L (x ,t ) H 2 * n Lemma 2.2 are used, 
and 



I 1 (x ,t ) = 2V2c 2 [ f \u\ 2 \v\ 2 dxdt + V2c f X ° ° \u \ 2 dx, 
I 2 (x ,t ) = 2V2c 2 f f \u\ 2 \v\ 2 dxdt + \f2c \v \ 2 dx. 

J JQtxn.tn) Jxn—tn 



>n(xo,to) Jxo-to 

Therefore these lead to the desired result by Lemma 2.1. The proof is complete. 
□ 

Now, to get the H 1 bound of the solution, we differentiate the equations 
(1.1) with respect to x, then 

f i(u xt + u xx ) = V u ,v,u,vGi(u,v) ■ (u x ,v x ,u x ,v x ), 
\ i(v x t -v xx ) = V u , v ,u,vG 2 {u,v) ■ (u x ,v x ,u x ,v x ), 

which lead to the following 

sup (\\u(t)\\ H i + \\v(t)\\ H i) < (\\u \\ H i + \\v \\ H i)exp(C%) < oo, 

0<t<T* 



7 



where the constant C depends only on L°° bounds of (u, v) given in Lemma 
2.3. Therefore, the local solution could be extended globally. Since C^°(R 1 ) is 
dense in with the above estimates, we get the following. 

Proposition 2.1 Suppose that (u ,v ) G and that ||mo||l 2 (.r 1 )+|| 1 'o||.l 2 (.r 1 ) 

is small enough. Then (1.1-1.2) has a unique global classical solution (u,v) G 
C , ([0,+oo);/J 1 (i? 1 ))nC 1 ([0,+oo);L 2 ( J R 1 ))nL oo ( J R 1 x [0,oo)). Moreover, there 
exists some positive constant C 2 such that 

sup (\u(x,t)\ 2 + \v(x,t)\ 2 ) < C 2 sup (\u (x)\ 2 + \v (x)\ 2 ), 

t>0,xeR} xeR 1 

[ b (\u(x,t)\ 2 + \v(x,t)\ 2 )dx < C 2 [ b (\u (x - t)\ 2 + \v (x + t)\ 2 )dx, 

J a J a 

for any a < b. 



3 Proof of the main result 

Let (uj t o, Vjfi) G C^°(R 1 ) be a sequence of smooth functions such that (tij,o, Vj,o) — > 
(uo,Vq) strongly in L 2 and such that sup^- ||(wj,o, ^,0)1^°° < 00. 

Denote (uj,Vj) be the if ^solutions of (1.1) with initial data (u^o, Vj,o)- 
Then by Proposition 2.1 we derive the following 

sup (| \uj(t)\\ L oo + \ \vj(t)\\ L °o) < 00. (3.1) 

t>0,j 

Now for any k,j 

(u k - Uj) t + (u k - Uj) x ) = iGi(uj,Vj) - iGi(u k ,v k ), 
(v k - Vj) t - (v k - Vj) x ) = iG 2 (uj, Vj) - iG 2 (u k , v k ), 

which together with (3.1) gives the following by the energy estimate method 

HK-^)WII! 2 + IIK-^)WII! 2 

<exp(O(l)t)(||K- % )(0)||l 2 + ||K-^)(0)||! 2 ), 

where the bounds of 0(1) depend only on the L°° bounds of the solutions 
(uj,Vj) (j = 1,2, • • •) for t > 0. This yields the convergence of the sequence 
{(«j,Uj)}£Lo to a function (u,v) in C([0, T); L 2 (R 1 )) for any T > 0. Then 
(u,v) is a global strong solution to (1.1-1.2). Moreover, by using Proposition 
2.1 and by repeating the argument of taking limit, we have 

[ b (\u(x,t)\ 2 + \v(x,t)\ 2 )dx < d [ b (\u (x - t)\ 2 + \v (x + t)\ 2 )dx 

J a J a 

rb+t rb-t 

= d( / \u (y)\ 2 dy+ / \v (y)\ 2 dy). 
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Then the decay of the local charge follows. 

In the same way, we can prove the following for any strong solutions (u', v') 
and (u",v"), 

1 1 - ^) (*) 1 1 + 1 ] - ^'0 (0 1 1 < exp (o(i)t) ( j | - ^'0 (o) ! | i 2 + 1 1 - ^) (o) 1 1 

for t > 0, where the bounds of 0(1) depend only on the L°° bounds of the 
solutions. This leads to the uniqueness of the strong solution. The proof of 
Theorem 1.1 is complete. 
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